The Stellar Opacity 


The mean absorption coefficient, K, is not a constant; it is de- 
pendent on frequency, and is therefore frequently written as «p. 
Inside a star, several different sources of opacity are important, 
and each has its own frequency dependence. To compute the ap- 
propriate “average” opacity, k recall that Fick’s law states that 
the flux at any point is related to the energy density U by 
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and that, when integrated over all energies, 
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But how is & related to k? To answer this question, we start by 
explicitly writing in the frequency dependence of U, i.e., 
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and compute the monochromatic flux distribution 
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The total flux, integrated over all frequencies is then 
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Comparing (6.1.1) and (3.1.4), we see that 
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The average opacity that you should use is therefore that which 
is harmonically weighted over the temperature derivative of the 
Planck curve. This is called the Rosseland mean opacity. 
Note that the function weights the high frequencies more than 
the low frequencies. To see this, we can take the temperature 
derivative of the Planck function 
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where x = hv /kT, and find its peak by taking the derivative. 
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Mathematically, this works out to x ~ 3.83. In other words, the 
high energy photons count much more than low energy photons. 


Although in real models, the Rosseland mean opacity is com- 
puted numerically, it is instructive to evaluate the form of & 
when ky xv”. From the definition of the Rosseland mean 
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Thus, if the opacity law is k x v`”, the temperature dependence 
of the Rosseland mean opacity will be k x T7”. This follows 
since the mean opacity must drop as the number of high energy 
photons increases. 
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ELECTRON SCATTERING 





When free electrons have thermal motions that are much less 
than their rest mass energy, i.e., Mec? > kT, or T < 5.9 x 
10° K, the cross-section for scattering of photons off electrons is 
frequency independent. The opacity per unit mass of material in 
this case is 
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where ce is the Thomson cross-section of the electron 
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and ne is the density of electrons 
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If the gas is fully ionized, i.e., if the ionization fractions f; ~ 1 (as 
it will be in stellar interiors), and if we approximate the atomic 
weight (A;) of a metal as twice its atomic number (Z;), then 
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where X, Y, and Z are the mass fractions of hydrogen, helium, 
and metals, respectively. Since X +Y + Z = 1, 
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and 


Ke = NACe(l + X)/2 = 0.2 (1 + X) cm? g! (6.1.7) 


Note again that as long as the ionization is complete (and matter 
is not degenerate), the result is independent of temperature and 
density. 


FREE-FREE ABSORPTION 


The energy emitted per unit volume by a single electron scattered 
by ions of charge Z (as a function of velocity, v) is 
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where e is the charge of the electron, and gg is the gaunt factor. 
To get the emissivity of thermal bremsstrahlung as a function of 
frequency, we integrate this over a Maxwellian distribution, noting 
that to get a photon with frequency v, the incident electron must 
have m,v7/2 > hv. Thus, Vmin = (2hv/m.)\/? and 
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Next, to go from emission to absorption, we can use the fact that 
in thermodynamic equilibrium, Kirchoff’s law must hold, hence 
the source function is S, = j,/a, = B,(T), where a,, is the 
absorption coefficient per electron (with units of inverse distance). 
Thus 
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or, upon substituting p = neue/Na and p=nypr/Na 
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The first term in the equation gives the total absorption coefficient; 
the exponential term at the end represents a negative contribution 
from stimulated emission. 





If we combine the constants, letting 
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To calculate the Rosseland mean free-free opacity, the inverse of 
this expression, weighted by the derivative of the Planck Function, 
must be integrated over all frequencies. This is tedious, but note 
that if we neglect the effects of stimulated emission, the frequency 
dependence of the opacity is y~°, which, through (6.1.5), implies 
a temperature dependence of T~?. This, plus the T712 depen- 
dence which is already in the equation, gives a total temperature 
dependence of T~7/2. Numerically, 
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Opacities of the form k x pT are called Kramers opacities. 


BOUND-FREE ABSORPTION 


For a hydrogenic atom, the cross section for absorption from state 
(n,l) to the continuum is 
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where the energy of the incident photon must be great enough to 
ionize the atom, t.e., 
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From this, it is clear that the opacity law for a single species of 
atom will have a number of absorption “edges” below which the 
absorption cross section will drop quickly. The total bound-free 
opacity, of course, will be the sum over all elements, Z, all ioniza- 
tion states, 7, and and all excitation states, n. In the hydrogenic 
ion approximation, the opacity is 
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This summation, along with the process of taking the Rosseland 
mean, will act to smear out many of the edges. 


The bound-free opacity must be computed numerically, but we can 
get a rough idea of how bound-free opacity behaves by making the 
assumption that in a star’s interior, all ions will be very highly 
ionized and thus will have at most one electron. The Saha equation 
states that the number of atoms in ionization state į is related to 
the number in ionization state t + 1 by 
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where u(T) are the partition functions and x; is the ionization 
energy. If the gas is highly ionized, Nj; ~ 1, ui4a = 1, and 
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We can then estimate the number of these atoms in excitation 
state n from the Boltzmann equation. For hydrogenic ions, the 
statistical weight of a level is n°, so if we let Xin = Xi + Xn, then 
the total number of absorbing ions is 
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The total absorption coefficient is therefore 
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and the opacity is 


an ; ZEK 
emaon So (pert) (Ta) TE 


where we have dumped the many constants into the variable C 
(and have ignored the stimulated emission term). 


Now examine the factor x;,,/kT. If for a particular n and Z, 
Xin < kT, then that term will contribute little to the opacity, 
since the species will have too few bound electrons. Similarly, if 
Xin > kT, there will not be many photons of the requisite energy 
to be absorbed, and the contribution to the opacity will again be 
small. Thus, only those terms with Xin ~ kT will contribute. 
This leaves kx x pT~'/2v—%, which is the same as free-free opac- 
ity. Bound-free absorption will therefore also follow the Kramers 
opacity law. 





A rough numerical estimate for the Rosseland mean for bound-free 
opacity is 
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Note that this is a factor of ~ 100 greater than the Rosseland 
mean for free-free opacity. 
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To first order, the cross-section for ionization of most species goes 
as v`’ for frequencies above the ionization energy. 


BOUND-BOUND ABSORPTION 


At extremely high temperatures, T > 10° K, the radiation field will 
consist mostly of ionizing photons, hence bound-bound absorption 
will be minimal (< 10%). However, at cooler temperatures the 
onset of bound-bound transitions in the UV and far UV can double 
the stellar opacity. 


Obviously, the computation of bound-bound opacity must be done 
numerically, since it involves the summation of millions of individ- 
ual absorption lines. It should be noted, however, that in the 
interior of a star, the pressure (and thermal) broadening terms 
will be substantial, so the net result will be continuum absorp- 
tion, rather than a series of absorption lines. Also, the expression 
for the absorption coefficient of hydrogenic bound-free transitions 
(6.1.12) is the same as that for bound-bound transitions. Thus, we 
can expect the same ~ v~? dependence and therefore a Kramer- 
like opacity law. The strength of the total absorption, however, 
will be critically dependent on temperature. 








H- OPACITY 


A free electron passing by a neutral hydrogen atom can produce 
a dipole moment in the atom. At that time, an extra electron can 
attach itself to the system; the result is an H~ ion. This ion is 
fragile — it has an ionization energy of only 0.754 eV, and a large 
cross section for absorption, either through bound-free absorption 
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In the envelopes of cool stars, H~ can be the dominant source of 


opacity. 


The contribution of H~ to the star’s opacity will depend on the 
number of H~ ions, which can be calculated from the Saha equa- 
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where upg- = 1, up, = 2, and y = 0.754 eV. If x is the fraction 
of ionized hydrogen (x ~ 0 for cool stars), and a, is the frequency 
dependency of H~ absorption, then the opacity due to H™ is 
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Note that the total opacity is directly proportional to the electron 
density. In extremely metal poor stars, the free electrons must 
come from ionized hydrogen, hence the opacity will rise with tem- 
perature, until the number of neutrals starts to decline. However, 
for normal metal rich stars, common, low ionization metals (such 
as Ca, Na, K, and Al) contribute to the supply of free electrons, 
thereby lessening this dependence. Consequently, H~ opacity is 
much stronger in metal rich stars than in metal poor stars, and is 
stronger in dwarfs than in giants (due to the pressure term). 














The determination of the frequency dependent cross-section of the 
H- ion involves the quantum mechanical computation of overlap- 
ping wave functions. It is not a power law: it increases linearly at 
short wavelengths until ~ 8500 A, then declines until the ioniza- 
tion threshold at ~ 16,500 A. After that, that opacity increases 
again due to free-free opacity. As a result, the Rosseland mean of 
H~ opacity is not a Kramers law opacity. Very roughly 
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Bound-free and free-free opacity from H7 at T = 6300°K. Ordinate: 
cross-section { x 104°] per neutral H atom and per unit electron 
pressure p, = AAT; abscissa: 4/1000 where 2 is in A, 


The combined bound-free and free-free opacity from the H~ ion 
at T = 6300 K. The y-axis plots the cross-section (x107~?° per 
neutral H atom per unit electron pressure, Pe = NnekT). From 
Mihalas, Stellar Atmospheres (2nd edition), pg. 103. 
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This figure illustrates the importance of the various opacities, as 
a function of temperature and density, for a normal (Population 
I) mix of abundances. The lines are labeled in units of electron 
scattering opacity, Ko = Ke = 0.2(1 + X). From Hayashi, Hoshi, 
& Sugimoto 1962, Suppl. Prog. Theo. Phys. Japan, Vol. 22. 
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The Los Alamos National Labs radiative opacities for X = 0.7, 
Y = 0.28, and Z = 0.02. The dashed line shows the half-ionization 


curve for pure hydrogen. 
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The Rosseland mean opacity (in cm? g~') as a function of density 
(in g cm~?) and temperature for X = 0.739, Y = 0.24, and Z = 
0.24. The range of values plotted is applicable for the outer regions 
of stars; the dots represent values in the solar model. 


Tabulated Opacities 


When computing real stellar models, simple power-law represen- 
tations of the Rosseland mean opacity are inadequate. Instead, 
one uses tables created either by The Opacity Project 


http://cdsweb.u-strasbg.fr/topbase/The0P.htm1 
or OPAL 
http://opalopacity.11n1.gov 


One chooses the table appropriate for the metallicity being consid- 
ered, and then interpolates in density and temperature (or related 
variables). Note that most of these tables assume a fixed metal 
distribution (i.e., the amount of given element, such as oxygen, 
scales with that of iron). However, a few do allow for independent 
changes in CNO, or a-process elements. 
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Opacity versus temperture for various measures of R = p/T;. The 
plots are for solar composition of 17 elements, including H, He, C, 
N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe, and Ni. From 
the Opacity project. 


The Eddington Luminosity 


If energy transport at the surface of a star is through radiation, 
then equation (3.1.4) holds, and the star’s luminosity is related to 
its temperature gradient by 
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Now assume that the dominate source of pressure is from radiation. 
(At the surface of a bright star, this will certainly be true.) For 
radiation pressure, 
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This is the Eddington luminosity. If the luminosity of the star is 
greater than this, the surface cannot be stable, and will be ejected 
from the star via radiation pressure. 


Note that equation (6.2.6), coupled with the mass-luminosity re- 
lation on the main sequence, implies that there is an upper limit 
to the mass of main sequence stars. If, for main sequence stars, 
Lr x M?-° and electron scattering is the dominant source of opac- 
ity at the surface, then 
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This equation implies that the maximum mass of a main-sequence 
star is Mmaz ~ 60Mo. Stars more massive than this must have 
strong stellar winds and will lose mass rapidly. 


